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Abstract 

In some options for circular polarization control at X-ray FELs, a helical radiator is 
placed a few ten meters distance behind the baseline undulator. If the microbunch 
structure induced in the baseline (planar) undulator can be preserved, intense 
coherent radiation is emitted in the helical radiator. The effects of betatron motion 
on the preservation of micro bunching in such in-line schemes should be accounting 
for. In this paper we present a comprehensive study of these effects. It is shown that 
one can work out an analytical expression for the debunching of an electron beam 
moving in a FODO lattice, strictly valid in the asymptote for a FODO cell much 
shorter than the betatron function. Further on, numerical studies can be used to 
demonstrate that the validity of such analytical expression goes beyond the above- 
mentioned asymptote, and can be used in much more a general context. Finally, a 
comparison with Genesis simulations is given. 



1 Introduction 

The LCLS baseline includes a planar undulator system, w^hich produces in- 
tense linearly polarized light in the wavelength range between 0.15 nm and 
1.5 nm ijlt]. Several schemes using helical undulators have been proposed for 
polarization control at the LCLS setup . The option presented in [4], 

exploits the microbunching of the planar undulator. After the baseline un- 
dulator, the electron beam is transported along a 40 m long straight line by 
FODO focusing system and subsequently passed through a helical radiator. 
If the microbunch structure of the bunch can be preserved, intense coherent 
radiation is emitted in the helical radiator. The driving idea of this proposal 
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is that the background linearly-polarized radiation from the baseline un- 
dulator is suppressed by spatial filtering. This operation consists in letting 
radiation and electron beam through horizontal and vertical slits upstream 
of the helical radiator, where the radiation spot size is about ten times larger 
than the electron bunch transverse size. The effect of betatron motion on the 
preservation of micro bunching in such scheme should be accounting for. 
In fact, the finite angular divergence of the electron beam, linked with the 
betatron function, yields a spread of the longitudinal velocity leading to mi- 
crobunching suppression. In 10] we estimated this factor, and concluded that 
the betatron motion should not constitute a serious problem in proposed 
scheme. 

In this paper we present a comprehensive study of the effect of betatron 
motion on the microbunch preservation. Our paper is based on the use 
results of llsj], where it was showed that in the limit for a small length of 
a FODO cell with respect to the betatron function value, the longitudinal 
velocity of an electron, averaged over a FODO cell, is constant through the 
focusing system. Based on this non-trivial statement, one can work out an 
analytical expression for the debunching of an electron beam moving in a 
FODO lattice, strictly valid in the asymptote for a short FODO cell. Further 
on, numerical studies can be used to demonstrate that the validity of such 
analytical expression goes beyond the above-mentioned asymptote, and can 
be used in much more a general context. 

The present work is organized as follows. In the following Section 2 we 
review the main result in [5*], and we report the expression for the average 
longitudinal velocity of an electron, which depends on the Courant-Snyder 
invariant of motion. Based on this result, the analytical expression for the 
debunching is calculated. In Section 3, the analytical asymptote is cross- 
checked with numerical calculations, and its validity is extended. A com- 
parison with results obtained with the code Genesis [|6fl is presented. Finally, 
in Section 4, we come to conclusions. 



2 Analytical study 



The starting point for our analysis is a consideration on the longitudinal 
velocity of an electron in a FODO lattice. As has been shown in ijij], when the 
length of the FODO cell LpoDO is much shorter than the betatron function, the 
longitudinal velocity of the electron, averaged over one FODO cell length, is 
constant. From a mathematical viewpoint, this result may be obtained with 
simple analytical calculations from Eq. (6) and Eq. (7) of reference [5]. We will 
begin to consider a 2D motion on the x(horizontal) - z(longitudinal) plane. 
The longitudinal velocity of a certain electron can be written as Uz(z) = v[l - 



3 



cos(0(z))], where v is the electron speed, and d{z) = x'{z) is the angle formed 
at each point of the electron trajectory with the longitudinal axis. When, as 
is the case for ultrarelativistic electrons moving along the z axis, x'(z) ^ 1, 
one can expand the trigonometric function and obtain v{z) — v - ux'(z)^/2, 
where x'{z) is given by Eq. (6) and Eq. (7) of reference When the length 
of FODO cell, Lfodo/ is much shorter then betatron function, the magnitude 
of the Twiss parameter a approaches unity, and one can re-write these two 
equations approximatively as: 



Here Ix is the particle Courant-Snyder invariant, while j6 and (p are the 
betatron function and the betatron phase respectively Averaging over 
one FODO cell length one obtains: 



which is proportional to the Courant-Snyder invariant I^, and is independent 
of z. This result looks at first glance surprising. In fact, a glance of a typical 
electron trajectory shows. Fig. 1 shows an overall oscillatory trajectory, and 
one would expect that the longitudinal velocity at z = 0, where the electron 
trajectory forms an angle with the z axis, should be smaller than, for example, 
that at z - 10 m, where the trajectory is almost parallel to the z axis. The 
physical explanation is that the longitudinal velocity is decreased at the 
position for at z ^ 10 m, by the presence of sharp oscillations on the scale 
of the FODO cell, which, as demonstrated above, lead to a constant average 
longitudinal velocity. 

We will now make use of Eq. (2) to study the effects of the betatron motion on 
the preservation of EEL microbunching. Consider an electron beam carrying 
an average current Iq. Let us superimpose an initial modulation at a given 
frequency co. The total current can be written as a function of the phase 
■ip = oj{z/ Vz-t), with z the longitudinal position, t the time, Vz the longitudinal 
velocity as 




(1) 




(2) 







Ii = 7o(l + fli cos ip) . 



(3) 
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Fig. 1. Trajectory of an electron within a FODO cell. Here LpoDO = 1-0 m and j3 ^ 10 
m. 

Similarly, given the unmodulated longitudinal particle density Hq = Io/{-eVz), 
{-e) being the electron charge, the longitudinal particle density ni{xjj) after 
the initial modulation is given by 



ni = no(l + fli cos i/') . (4) 

Here fli describes the initial bunching. The relation between Ui and the 
bunching factor bi =< exp{iip) >, which can be often found in literature is 
given by 



Nptc 27T 



^1 = 2-1 ^^P^^^^'^ "in J ^^~n^ exp(zi/') = — , 

k=l r\ 



(5) 



where Nptc is the total number of particles within a wavelength A = InVz/co, 
i.e. Nptc - WqA, and ipk is the phase of each particle. 

The ratio L^odoIP> with j3 the average betatron function, is the first rele- 
vant parameter of our problem. We will assume LpoDo/jS ^ 1, so that Eq. 
(2) can be used for the average longitudinal velocity of an electron. The 
phase difference of an electron with Courant-Snyder invariant Ix with re- 
spect to one moving on-axis, after a given distance L can thus be written as 
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dip = AvzL/{2cA) ^ -IxL/{2^A). Since the rms value for Ix is the geometrical 
emittance e one obtains, parametrically, dip is of order eL/(j6A), which is 
the second and last relevant parameter of our problem. No debunching is 
expected for £L/(jSA) ^ 1. 

Under the accepted approximation Lfodo/iS <^ 1 it is possible to derive 
an analytical expression for the evolution of the bunching factor along the 
FODO lattice. In fact, the influence of the betatron motion alone can be 
modelled by substituting the phase xpk of each individual electron with 
i/^fc + 6ip]^{dx) where, as described above. 



dipk = 



(6) 



The bunching factor &2 after the propagation along the FODO cell can there- 
fore be written as an average over the distribution of Ix, which we call f{Ix), 
as 



Nptc 



&2 =< exp[z(i/'fc + dipk)] >= 2^ exp[i{ipk + 6xpk)] 

k=l 

2n oo 

= J dlxf{lx)exY>{i{ip + bip)]. 



Using 



(7) 



fi^x) = 2^ exp 



k 

2e 



(8) 



remembering Eq. (4) and Eq. (5) one obtains from Eq. (7) the following 
expression for the ratio C between final and initial bunching: 



b2 1 r,. 

C, = — ^ — ^ — dixexp 
ai bi 2e J 



- Ix' 




ilxL 


. 2£. 


exp 


2j3A 



(9) 



The integral can be calculated analytically yielding the final result 



(10) 



Note that C is not a real number. The physical interpretation of \Q = [1 + 
£^L^/(A^j8^)]"^''^ is that of the evolution of the amplitude of the bunching. 
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Fig. 2. Modulus (left plot) and phase (right plot) of C as a function of the parameter 

The physical interpretation of Arg(C) is that of the evolution of the phase of 
the bunching. In other words, the bunching evolves in modulus and phase 
along the FODO lattice. Also, notice that C is independent of the initial 
definition of the bunching, ui or bi. In both cases, the final bunching can be 
found by multiplying the initial bunching by C- Modulus and phase of C are 
plotted in Fig. 2 as a function of the parameter eL/ (j3A). 

Finally, we note that all considerations have been made for a 2D motion on 
the x-z plane. However, generalizing to a 3D motion is trivial. One needs 
to account the divergence y', giving an extra contribution to the phase 
formally identical as that for the horizontal direction, once the horizontal 
beta function is substituted with the vertical one. One obtains that C = CxCy, 
with Cx,y formally identical to Eq. (10). 



3 Numerical study 

In order to study the influence of the betatron beating on the microbunch 
suppression and to have an idea about the accuracy of our analytical asymp- 
totic results, we simulated the evolution of the bunching numerically. In par- 
ticular, we considered a periodic lattice composed of drift, focusing element 
in the thin-lens approximation, drift, def ocusing element. For simplicity, we 
considered the motion in the x-z direction, so that Eq. (10) could be used to 
compare with debunchtng calculated numerically. In the numerical calcula- 
tions, the magnetic structure is defined, in terms of quadrupole strength and 
length of the drifts. Then, the Twiss parameters at the beginning of the setup 
are calculated, and used to generate the horizontal phase-space distribution 
of electrons. Each particle is tracked through the setup in a linear matrix 
approach, and the trajectory is calculated. Knowing the trajectory Xk{z) for 
each electron it is then straightforward to calculate the curvilinear distance 
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Fig. 3. Comparison between |C| calculated analytically as a function of the emittance 
from Eq. (10) (solid line) and calculated numerically (black circles) as described in 
this Section. Here j3 ^ 10m, and LpoDO = liri- The modulation wavelength is A = 1.5 
nm, while the drift length is L = 40 m. The number of particles is varied from 10^ to 
5-10^ to achieve a good accuracy. As an exemplification, the white circle corresponds 
to 10^ particles, whereas the black circle at the same emittance value corresponds 
to 5 • 10^ particles. 

traveled. Comparison with the length of the setup allows one to recover the 
phase difference 6ipk and to calculate C =< exp{i6\l^k) >• 

We began our numerical investigations setting the average betatron function 
value jS - 10 m, and the length of the cell LpoDO = liri. This choice for the 
length of the FODO cell allows to consider the asymptote LpoDo/iS <c 1 
nearly satisfied. The debunching module \Q was calculated for different 
emittance value. A comparison with |CI calculated analytically from Eq. (10) 
is shown in Fig. 3. The number of particles is varied from 10"* to 5 • 10^ to 
achieve a good accuracy, and is not constant for the calculated points. For 
the sake of exemplification, the white circle in Fig. 3 corresponds to 10^ 
particles, whereas the black circle at the same emittance value corresponds 
to 5 • 10^ particles. The good agreement between numerical and analytical 
results was to be expected on the basis of the asymptote L^odo/^ 1- 

An interesting result can be achieved by fixing the emittance (in our case 
£ = 1.188 • 10"^° m), and changing j8, keeping the other quantities as in 
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Fig. 4. Comparison between \C,\ calculated analytically as a function of the average 
betatron function from Eq. (10) (solid line) and calculated numerically (black circles) 
as described in this Section. Here e - 1.188 • 10"^*^ m, and LpoDO = Im- The 
modulation wavelength is A = 1.5 nm, while the drift length is L = 40 m. 

the previous example. This means that still LpoDO = 1 m. In this way one 
can sweep through different values of the parameter LpoDo/jS- In particular, 
we calculated numerically the value of |C| for 1/40 < LpoDo/jS < 1/2. The 
result of the comparison with the asymptotic for L^odqI^ «; 1 in Eq. (10) 
is shown in Fig. 4. A very good agreement can be seen even for values 
i-FODo/jS not much smaller than unity. Since the average beta function is 
related to the betatron phase (p and to the length of the FODO cell by 
jS = LFODo[cot((/)/2) + 2/3tan((/)/2)], the maximuum value achievable for 
^FODo/jS turns out to be limited to about 0.61 for a working focusing system. 
From this numerical analysis it follows that Eq. (10) can be used not only in 
the asymptotic case for Lfodo/iS ^ 1 but, with the good accuracy given in 
Fig. 4, practically in all cases. This is the main result of this work. 

Finally, we present a comparison between |C| obtained from Eq. (10), cal- 
culated numerically as above, and derived using the EEL code Genesis [6]. 
We assumed a cell length of 3.84 m, so that we set Lpooo = 3-84m, and a 
drift distance equivalent to 10 cells, i.e. L = 38.4 m. We set both horizontal 
and vertical average betatron function to 10.1 m. In order to simulate the 
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Fig. 5. Evolution of the rms horizontal and vertical beam size as a function of the 
distance along the setup calculated through Genesis at a normalized emittance 
£ „ = 4 ■ 10'' m and at an electron beam energy of 4.3 GeV (y = 8416). 

focusing system in Genesis, without the influence of the undulators we set 
the electron beam current to zero, we switched off the undulator focusing, 
and we prepared a Genesis particle file with a given density bunching at 
the modulation wavelength A = 1.5 nm, so that it was matched with the 
FODO beam transport line. The electron beam energy was set to 4.3 GeV. 
All particles in the particle file were set with the same energy: as a result 
effects of the momentum compaction factor were excluded. The beam was 
propagated through the setup. The evolution of the rms horizontal and ver- 
tical size as a function of the distance along the setup is shown in Fig. 5. At 
the end of the setup, the final particle beam was extracted, allowing for a 
comparison of the final bunching with respect to the initial bunching. The 
debunching as a function of the geometrical emittance is presented in Fig. 
6. 



4 Conclusions 



In this paper we derived an analytical expression for the debunching of 
an modulated electron beam through a FODO focusing structure. The ex- 
pression is very simple, and can be practically used for any value in the 
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Fig. 6. Comparison between \C,\ calculated analytically as a function of the geomet- 
rical emittance from Eq. (10) (solid line), |CI calculated numerically (black circles) 
as described in this Section and \Q calculated through Genesis. 

parameter space. 
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